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Abstract 

Spectroscopic analyses with the intention of the interpretation of the UV-spectra 
of the brightest stars as individuals - supernovae - or as components of star- 
forming regions - massive O stars - provide a powerful tool with great astro- 
physical potential for the determination of extragalactic distances and of the 
chemical composition of star-forming galaxies even at high redshifts. 

The perspectives of already initiated work with the new generation of tools 
for quantitative UV-spectroscopy of Hot Stars that have been developed during 
the last two decades are presented and the status of the continuing effort to 
construct corresponding models for Hot Star atmospheres is reviewed. 

Since the physics of the atmospheres of Hot Stars are strongly affected by 
velocity expansion dominating the spectra at all wavelength ranges, hydrody- 
namic model atmospheres for 0-type stars and explosion models for Supernovae 
of Type la are necessary as basis for the synthesis and analysis of the spectra. 
It is shown that stellar parameters, abundances, and stellar wind properties can 
be determined by the methods of spectral diagnostics already developed. Ad- 
ditionally, it will be demonstrated that models and synthetic spectra of Type la 
Supernovae of required quality are already available. These will make it pos- 
sible to tackle the question of whether Supernovae la are standard candles in 
a cosmological sense, confirming or disproving that the current SN-luminosity 
distances indicate accelerated expansion of the universe. 

In detail we discuss applications of the diagnostic techniques by example of 
two of the most luminous O supergiants in the Galaxy and a standard Supernova 
of Type la. Furthermore, it is demonstrated that the spectral energy distributions 
provided by state-of-the-art models of massive O stars lead to considerably bet- 
ter agreement with observations if used for the analysis ofHll regions. Thus, an 
excellent way of determining extragalactic abundances and population histories 
is offered. 

Moreover, the importance of Hot Stars in a broad astrophysical context will 
be discussed. As they dominate the physical conditions of their local environ- 
ments and the life cycle of gas and dust of their host galaxies, special emphasis 
will be given to the corresponding diagnostic perspectives. Beyond that, the 
relevance of Hot Stars to cosmological issues will be considered. 
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1 Introduction 

It is well known that Hot Stars are not a single group of objects but comprise sub- 
groups of objects in different parts of the HR diagram and at different evolution- 
ary stages. The most important sub-groups are massive O/B stars, Central Stars of 
Planetary Nebulae, and Supemovae of Type la and II. All these sub-groups have in 
common that they are characterized by high radiation energy densities and expand- 
ing atmospheres, and due to this the state of the outermost parts of these objects is 
characterized by non-equilibrium thermodynamics. In order to cover the best-known 
fundamental stages of the evolution of Hot Stars in sufficient depth this review will 
be restricted to the discussion of O stars and Supernovae of Type la; we will not dis- 
cuss objects like Wolf-Rayet stars. Luminous Blue Variables, Be-stars, Supemovae 
of Type 11, and others. Furthermore, Hot Stars play an important role in a broad 
astrophysical context. This impUes that a complete review covering all aspects in 
theory and observation is not only impossible, but also beyond the scope of this re- 
view. Thus, this review will focus on a special part of the overall topic with the intent 
to concentrate on just one subject; the subject chosen is UV spectral diagnostics. It 
will be shown, however, that this subject has important implications for astrophysical 
topics which are presently regarded as being "trendy". 

But first of all we have to clarify what UV spectral diagnostics means. This is best 
illustrated by the really old-fashioned (1977, Morton and Underbill) UV spectrum of 
one of the brightest massive O stars, the 04I(f) supergiant C, Puppis. As can be 
seen in Figure 1, expanding atmospheres have a pronounced effect on the emergent 
spectra of hot stars - especially in the UV-part. The signatures of outflow are clearly 
recognized by the blue-shifted absorption and red-shifted emission in the form of the 
weU-known P Cygni profiles. It is quite obvious that these kind of spectra contain 
information not only about stellar and wind parameters, but also about abundances. 
Thus, in principle, all fundamental parameters of a hot star can be deduced from a 
comparison of observed and synthetic spectra. 

Although spectra of the quahty shown in Figure 1 have been available for more 
than 25 years, most of the work done to date has concentrated on qualitative results 
and arguments. In view of the effort put into the development of modern - state-of- 
the-art - instruments, it is certainly not sufficient to restrict the analyses to simple 
line identifications and qualitative estimates of the physical properties. The primary 
objective must be to extract the complete physical stellar information from these 
spectra. Such diagnostic issues principally have been made possible by the superb 
quality and spectral resolution of the spectra available. 

For this objective, the key is to produce realistic synthetic UV spectra for Hot 
Stars. Such powerful tools, however, are still in development and not yet widely 
used. They offer the opportunity to determine the stellar parameters, the abundances 
of the light elements - He, C, N, O, Si - and of the heavy elements like Fe and Ni, 
quantitatively, as is indicated not only by lUE, but also more recent HST, ORFEUS, 
and FUSE observations of hot stars in the Galaxy and Local Group galaxies. All 
these observations show that the spectra in the UV spectral range are dominated by 
a dense forest of slightly wind-affected pseudo-photospheric metal absorption lines 
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Figure 1 : Merged spectrum of Copernicus and lUE UV tiigh-resolution observations of the 
04 1(f) supergiant C Puppis (900-1500 A: Morton and Underbill 1977; 1500-1800 A: Walbom 
etal. 1985). The most important wind Unes of the Ught elements are identified and marked. 
Also marked are the large number of wind-contaminated lines of the iron group elements (e.g., 
Fe V) which are especially present between 1250 and 1500 A. (Figure from Pauldrach etal. 
1994b). 



overlaid by broad P Cygni line profiles of strong, mostly resonance lines, formed 
in different parts of the expanding atmosphere. Thus, the obvious objective is to 
investigate the importance of these lines with respect to the structure of the expanding 
atmospheres that are characterized by the strength and the velocity of the outflow, 
and through which the shape of the spectral lines is mainly determined. Over the last 
30 years it turned out that the achievement of this objective, which will also be the 
primary aim of this paper to review, remains a difficult task. 

Before we discuss in detail the status quo of the diagnostic tool required (Sec- 
tion 6), we will first examine whether such a tool has been made obsolete by the 
general development in astrophysics or whether it is stiU relevant to current astro- 
nomical research. 

For this purpose we first discuss the diagnostic perspectives of galaxies with 
pronounced current star formation. Due to the impact of massive stars on their en- 
vironment the physics underlying the spectral appearance of starburst galaxies are 
rooted in the atmospheric expansion of massive O stars which dominate the UV 
wavelength range in star-forming galaxies. Therefore, the UV-spectral features of 
massive O stars can be used as tracers of age and chemical composition of starburst 
galaxies even at high redshift (Section 2). 
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With respect to the present cosmological question of the reionization of the uni- 
verse - which appeared to have happened at a redshift of about 2: ~ 6 - the ionization 
efficiency of a top-heavy Initial Mass Function for the first generations of stars is dis- 
cussed (Section 3). 

Starting from the impact of massive stars on their environment it is demonstrated 
that the spectral energy distributions provided by state-of-the-art models of massive 
O stars lead to considerable improvements if used for the analysis of H 11 regions. 
Thus, the corresponding methods for determining extragalactic abundances and pop- 
ulation histories are promising (Section 4). 

Regarding diagnostic issues, the role of Supemovae of Type la as distance indi- 
cators is discussed. The context of this discussion concerns the current and rather 
surprising surprising result that distant SNe la appear fainter than standard candles 
in an empty Friedmann model of the universe (Section 5). 

Finally, we discuss applications of the diagnostic techniques by example of two 
of the most luminous O supergiants in the Galaxy; additionally, basic steps towards 
realistic synthetic spectra for Supemovae of Type la are presented (Section 7). 

2 The impact of massive stars on their Environment - 
UV Spectral Analysis of Starburst Galaxies 

The impact of massive stars on their environment in the present phase of the universe 
is of major importance for the evolution of most galaxies. Although rare by number, 
massive stars dominate the life cycle of gas and dust in star forming regions and 
are responsible for the chemical enrichment of the ISM, which in turn has a signifi- 
cant impact on the chemical evolution of the host galaxy. This is mainly due to the 
short lifetimes of massive stars, which favours the recycUng of heavy elements in an 
extremely efficient way. Furthermore, the large amount of momentum and energy 
input of these objects into the ISM controls the dynamical evolution of the ISM. 
This takes place in an extreme way, because massive stars mostly group in young 
clusters, producing void regions around themselves and wind- and supernova-blown 
superbubbles around the clusters. These superbubbles are ideal places for further 
star formation, as numerous Hubble Space Telescope images show. Investigation 
of these superbubbles will finally yield the required information to understand the 
various processes leading to continuous star formation regions (cf. Oey and Massey 
1995). The creation of superbubbles is also responsible for the phenomenon of galac- 
tic energetic outflows observed in starbursts (Kunth et al. 1998) and starburst galaxies 
even at high redshift (Pettini etal. 1998). It is thus not surprising that spectroscopic 
studies of galaxies with pronounced current star formation reveal the specific spec- 
tral signatures of massive stars, demonstrating in this way that the underlying physics 
for the spectral appearance of starburst galaxies is not only rooted in the atmospheric 
expansion of massive O stars, but also dominated by these objects (cf. Figure 2 from 
Steidel etal. 1996, for star-forming galaxies at high redshift see also Pettini etal. 
2000, and for UV Une spectra of local star-forming galaxies see Conti etal. 1996; 
note that the similarity of the spectra at none/low and high redshifts suggests a simi- 
lar stellar content). 




Figure 2: UV spectrum of a z > 3 galaxy (upper part). For comparison, a recent HST 
spectrum of the central starburst region in the Wolf-Rayet galaxy NGC 4214 is also shown 
(lower part). Note the characteristic P-Cygni lines, especially of C IVand Si IV, pointing to the 
dominating influence of massive O stars. Figure from Steidel etal. 1996. 



The characteristic P-Cygni lines observed as broad stellar wind lines, especially 
those of the resonance lines of C IV and Si IV, integrated over the stellar populations 
in the spectra of starbursting galaxies, allow quantitative spectroscopic studies of 
the most luminous stellar objects in distant galaxies even at high redshift. Thus, in 
principle, we are able to obtain important quantitative information about the host 
galaxies of these objects, but diagnostic issues of these spectra require among other 
things synthetic UV spectra of 0-type stars as input for the population synthesis 
calculations needed for a comparison with the observed integrated spectra. 

The potential of these spectra for astrophysical diagnostics can nevertheless be 
investigated in a first step by using observed UV spectra of nearby O-type stars as 
input for the corresponding population synthesis calculations instead. In the frame of 
this method stars are simulated to form according to a specified star-formation his- 
tory and initial mass function and then follow predefined tracks in the HR-diagram. 
The integrated spectra are then built up from a library of observed UV spectra of 
hot stars in the Galaxy and the Magellanic Clouds. The output of this procedure are 
semi-empirical UV spectra between 1000 and 1800 A at 0.1 to 0.7 A resolution for 
populations of arbitrary age, star-formation histories, and initial mass function. The 
computational technique of this method is described comprehensively in the litera- 
ture and we refer the reader to one of the latest papers of a series (Leitherer etal. 
2001). 

As an example of the analyses performed in this way a comparison of the average 
spectrum of 8 clusters in NGC 5253 to synthetic models at solar (top) and 1/4 solar 
(bottom) metalUcity is shown in Figure 3 (Leitherer etal. 2001). The Figure shows 
clearly that a representative value of the overall metallicity of this starburst galaxy 
can be determined, since the model spectrum in the lower part fits the observation 
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Figure 3: Average spectrum of 8 clusters in NGC 5253 compared to that of a population 
model at solar (top) and 1/4 solar metallicity (bottom). It is clearly shown that a representa- 
tive value of the overall metallicity of this starburst galaxy can be determined. Figure from 
Leithereretal., 2001. 

almost perfectly. This result is especially convincing as both models (dashed lines) 
are based on the same parameters, except for the metallicity. A standard Salpeter 
IMF between 1 and 100 Mq was used and the starburst has been assumed to last 
6 Myr, with stars forming continuously during this time. The worse fit to the obser- 
vations produced by the solar metallicity model spectrum - particularly discrepant 
are the blue absorption wings in Si iv and C iv which are too strong in the models - 
could be improved by reducing the number of the most massive stars with a steeper 
IMF, but at the cost of the fit quality in the emission components. As an important 
result the ratio of the absorption to the emission strengths is therefore a sensitive 
indicator for the metalUcity. 

Moreover, the strong sensitivity of the emission parts of the P-Cygni hues of 
N V, Si IV, and C iv on the evolutionary stage of the O stars makes these spectra 
quite suitable as age tracers, which is shown in Figure 4, where the time evolution 
of the integrated spectrum following an instantaneous starburst is presented. The 
line profiles gradually strengthen from a main-sequence-dominated population at 
0-1 Myr to a population with luminous O supergiants at 3-5 Myr; and the lines 
weaken again at later age due to termination of the supergiant phase - this behavior 
is especially visible in the shape of the C iv A1550 resonance Une. 

The conclusion from these examples is that this kind of analysis is very promis- 
ing, but relies on observed UV spectra, which are just available for a small number 
of metallicity values, namely those of the Galaxy and the Magellanic Clouds. In or- 
der to make progress in the direction outUned before, realistic synthetic UV Spectra 
of 0-type stars are needed. This is in particular the case for high-redshift galaxies 
(which are observable spectroscopically when the flux is ampUfied by gravitational 
lensing through foreground galaxy clusters) since in these cases the expected metal- 
Ucities of starbursting galaxies in the early universe (cf. Pettini et al. 2000) are most 
probably different from local ones. 
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Figure 4: Time evolution of the integrated spectrum for tlie 10 Myr following an instanta- 
neous starburst. The strong sensitivity of the emission parts of the P-Cygni lines of N V, Si IV, 
and C IV on the evolutionary stage of the O stars makes these spectra quite suitable as age 
tracers. Figure from Leitherer etal. 2001. 

3 First generations of Stars - Ionization Efficiency 
of a Top-Heavy Initial Mass Function 

Apart from the short evolutionary timescale of massive stars it is obviously the metal- 
Ucity, and, connected to this, the steepness of the Initial Mass Function, which is 
responsible for the rarity of these objects in the present phase of the universe. It has 
to be the metalUcity, because very recently strong evidence has been found that the 
primordial IMF has favored massive stars with masses > lO^M© (cf. Bronnm etal. 
1999). Thus, in the early universe, when only primordial elements were left over 
from the Big Bang, nature obviously preferred to form massive stars. This predic- 
tion is based on the missing metallicity which leads to a characteristic scale for the 
density and temperature of the primordial gas, which in turn leads to a characteristic 
Jeans mass of Mj ~ WMq (Larson 1998). 

Finally, due to these physical conditions the Initial Mass Function becomes top- 
heavy and therefore deviates significantly from the standard Salpeter power-law (see, 
for instance, Bromm etal. 1999, 2001). Such an early population of very massive 
stars at very low metalUcities (i. e.. Population III stars) which have already been the- 
oretically investigated by El Eid et al. (1983), recently turned out to be also relevant 
to cosmological issues, the most important being the cosmological question of when 
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Figure 5: The normalized spectral energy distribution in the continuum of Population III 
stars - mass range 100-1000 Mq at Z — 0. Note that the spectra attain an universal form for 
stellar masses M > 300 Mq Figure from Bromm etal. 2001. 

and how the cosmic "dark ages" ended (Loeb 1998). We know that the dark ages 
ended because the absence of a Gunn-Peterson trough (Gunn and Peterson 1965) in 
the spectra of high-redshift quasars implies that the universe was reionized again at 
a redshift of z > 5.8 (Fan etal. 2000). For a long time, Carr etal. (1984) suspected 
that the first generations of stars have been relevant to control this process. Thus, 
Population HI stars could contribute significantly to the ionization history of the 
intergalactic medium (IGM), but the contribution of the first generation of stars to 
the ionization history of the IGM depends crucially on their initial mass function. 
With regard to the first generations of stars the ionization efficiency of a top-heavy 
Initial Mass Function will, therefore, have to be investigated. 

It is the enormous amount of UV and EUV radiation of these very massive stars 
which could easily change the status of the cold and dark universe at that time to 
become reionized again. This is indicated in Figure 5, which also shows that the total 
spectral luminosity depends solely on the total amount of mass, if the mass of the 
most massive stars exceeds 300 Mq (cf. Bromm etal. 2001). Due to this top-heavy 
Initial Mass Function the total spectral energy distribution deviates significantly from 
that obtained with the standard Salpeter power-law, as is shown in Figure 6, and the 
flux obtained can contribute the decisive part to the unexplained deficit of ionizing 
photons required for the reionization of the universe (cf. Bromm et al. 2001). 

However, in order to be able to make quantitative predictions about the influence 
of this extremely metal poor population of very massive stars on their galactic and 
intergalactic environment one primarily needs observations that can be compared to 
the predicted flux spectrum. 
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Figure 6: Predicted flux from a Population ni star cluster at 2 = 10. A flat universe with 
= 0.7 is assumed. The cutoff below Aobs = 1216 A is due to complete Gunn-Peterson 
absorption, and the observable flux is larger by an order of magnitude for the case of the 
top-heavy IMF when compared to the case of the standard Salpeter power-law. Figure from 
Brommetal. 2001. 

Future observations with the Next Generation Space Telescope of distant stel- 
lar populations at high redshifts will principally give us the opportunity to deduce 
the primordial IMF from these comparisons. Kudritzki (2002) has recently shown 
that this is generally possible, by calculating state-of-the-art UV spectra for massive 
O stars in a metallicity range of 1 . . . 10~^ Zq. From an inspection of his spec- 
tra he concluded that significant line features are stiU detectable even at very low 
metallicities; thus, there will be diagnostic information available to deduce physical 
properties from starbursting regions at high redshifts that will eventually be observed 
with the Next Generation Space Telescope in the infrared spectral region. 

As a second requirement we need to determine the physical properties of Popu- 
lation III stars during their evolution. A key issue in this regard is to obtain realistic 
spectral energy distributions calculated for metallicities different from zero for the 
most massive objects, since the assumption of a metalhcity of Z = OAs certainly 
only correct for the very first generation of Population III stars. 

4 Theoretical Ionizing Fluxes of O Stars 

Although less spectacular, we will now investigate the impact of massive stars on 
their enviroimient in a more direct maimer. Apart from the chemical enrichment of 
the ISM, the large amount of momentum and energy input into the ambient inter- 
stellar medium of these objects is primarily of importance. Especially the radiative 
energy input shortward of the Lyman edge, which ionizes and heats the Gaseous 
Nebulae surrounding massive Hot Stars, offers the possibiUty to analyze the influ- 
ence of the EUV radiation of the photoionizing stars on the ionization structure of 
these excited H 11 regions. 




10 



A. W.A. Pauldrach 




. . I I I I I I I 

0.01 0.10 1.00 10.00 

[Nelll]15.5/[Nell]12.8 



Figure 7: [Ann] 8.99 /im/[Arii] 6.99 /xm versus [Neiii] 15.6 /um/[Neii] 12.8 diag- 
nostic diagram of observed and predicted nebular excitation of H II regions. Boldface numbers 
indicate models which are designated by their effective temperature in 10^ K. Triangles are 
NLTE models and squares are LTE models; note that the NLTE models represent an impres- 
sive improvement in reproducing the [Ne III] emission. A representative error bar for the data 
is shown in the lower right comer Figure from Giveon et al. 2002. 

The primary objective of such investigations are studies of theoretical models 
of starburst regions, which for instance are used to determine the energy source 
in ultra-luminous infrared galaxies - ULIRGs - (cf. Lutz etal. 1996; Genzel etal. 
1998). The interpretation of the corresponding extra-galactic observations obviously 
requires understanding the properties of the spectral energy distributions (SEDs) of 
massive stars and stellar clusters. Thus, the quality of the SEDs has to be probed in 
a first step by means of investigations of Galactic H ii regions. 

4.1 The Excitation of Galactic H ii Regions 

Giveon etal. (2002) recently presented a comparison of observed [Nelll] 15.6 fim/ 
[Ne II] 12.8 and [Ar III] 8.99 /im/ [Ar 11] 6.99 /im excitation ratios, obtained for 
a sample of 112 Galactic H II regions and 37 nearby extragalactic H II regions in the 
LMC, SMC, and M33 observed with ISO-SWS, with the corresponding results of 
theoretical nebular models. 

The authors have chosen infrared fine-structure emission lines for their inves- 
tigation because these lines do not suffer much from dust extinction and the low 
energies of the associated levels make these lines quite insensitive to the nebular 
election temperature. JVloreover, the relative strengths of the fine-structure emission 
lines chosen are ideal for constraining the shape of the theoretical ionizing fluxes. 
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Figure 8: Comparison between the LTE (static atmosphere from Kurucz 1992) and NLTE 
(expanding atmosphere from Pauldrach etal. 2001) spectral energy distributions of a dwarf 
with an effective temperature of 35 000 K. The ionization edges of the relevant ions are in- 
dicated by vertical dashed lines. It is apparent that the difference in excitation produced by 
the LTE and NLTE atmospheres will be more pronounced in the Ne"'""'"/Ne"'" ratio than in the 
Ar+"'"/Ar"'" ratio. Figure from Giveon etal. 2002. 

since the [Nelll] 15.6 /im/[Nell] 12.8 /im line ratio depends on photons emitted 
at > 3 Rydbergs, while [Arm] 8.99 /zm/ [Arii] 6.99 jim is sensitive to the region 
> 2 Rydbergs. Thus, these line ratios are extremely useful probes of the physical 
properties of H II regions and their associated ionizing sources especially with regard 
to the ionizing spectral energy distributions. 

The complete set of observed emission line ratios of Ne and Ar is shown in 
Figure 7 together with the corresponding results of nebular model computations for 
which solar metallicity and a gas density of 800 cm~^ has been assumed, and which 
are based on LTE (static atmosphere from Kurucz, 1992) as well as NLTE (expanding 
atmosphere from Pauldrach etal., 2001) spectral energy distributions. 

The diagnostic diagram clearly shows that the predicted nebular excitation in- 
creases with increasing effective temperature of the photoionizing star, and it also 
shows that the high excitation [Ne III] emission observed in H II regions is by far not 
reproduced by nebular calculations which make use of the ionizing fluxes of LTE 
models - flie line ratios are under-predicted by factors larger than 10. This result is 
clearly an example of the well-known Ne iii problem (cf. Baldwin etal. 1991; Rubin 
etal 1991; Simpson etal. 1995). 

It is quite evident, however, that the quality of the diagnostics depends primarily 
on the quality of the input, i. e., the spectral energy distributions. It is therefore a 
significant step forward that the ionizing fluxes of the NLTE models used represent 
an impressive improvement to the observed excitation correlation - note that the fit 
shown in Figure 7 may actually be even better, since a lot of the scatter is due to un- 
derestimated extinction corrections for the [Ar III] 8.99 line. The NLTE sequence 
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also indicates that for most of the considered H II regions the effective temperatures 
of the exciting stars lie in the range of 35 000 to 45 000 K. Most importantly, the 
NLTE models can readily account for the presence of high excitation [Ne III] emis- 
sion Unes in nebular spectra. This result resolves the long-standing Ne iii problem 
and supports the conclusion of Sellmaier etal. (1996), who found, for the first time, 
on the basis of a less comprehensive sample of H II regions and models that this 
problem has been the failure of LTE photoionization simulations primarily due to a 
significant under-prediction of Lyman photons above 40 eV. This issue is solved by 
making use of NLTE model atmospheres. 

The reason for this improvement are the spectral shapes of the NLTE fluxes short- 
ward of the Ar II and Ne II ionization thresholds which are obviously somewhat more 
realistic in the NLTE case. This is illustrated in Figure 8 where the spectral energy 
distributions of an LTE and an NLTE model are compared by example of a dwarf 
with an effective temperature of 35 000 K. As is shown, the SEDs harden in the 
NLTE case, meaning that the LTE model produces much less flux above the Ne+ 
40.96 eV threshold than the NLTE model. This is the behavior that is essentially 
represented in the excitation diagram. 



4.2 Spectral Energy Distributions 
of Time-Evolving Stellar Clusters 

The improvement obtained for the diagnostic diagram of the analysis of Galactic 
H II regions in Section 4. 1 also has important implications for determining extra- 
galactic abundances and population histories of starburst galaxies. 

With regard to this, Thornley et al. (2000) carried out detailed starburst modelling 
of the [Neiii] 15.6 /xm/ [Neii] 12.8 /zm ratio of Hii regions ionized by clusters of 
stars. As was shown above, the hottest stars in such models are responsible for 
producing large nebular [Nelll] 15.6 /im/|Nelll 12.8 fim ratios. Hence, the low 
ratios actually observed led to the conclusion that the relative number of hot stars is 
small due to aging of the starburst systems. Thus, the solution of the Ne iii problem 
has important consequences for the interpretation of these extragalactic fine-structure 
line observations, since the conclusion that due to the low [Neiii] 15.6 /xm/ [Neii] 
12.8 /um ratios obtained, the hottest stars as dominant contributors to the ionization 
of the starburst galaxies have to be removed, has obviously to be proven. 

In order to tackle this challenge, Sternberg etal. (2002) computed spectral en- 
ergy distributions of time-evolving stellar clusters, on the basis of a large grid of 
calculated NLTE spectral energy distributions of O and early B-type stars (being the 
major contributors to the Lyman and He i fluxes) in the hot, luminous part (Teff > 
25 000 K) of the HR diagram. From this grid, models of stars following evolution- 
ary tracks are suitably interpolated. The SEDs used rely on recent improvements of 
modelUng expanding NLTE atmospheres of Hot Stars (Pauldrach etal. 2001). As 
an example of the individual models used. Figure 9 shows the calculated spectral 
energy distribution of a typical O star compared to the corresponding result of an 
NLTE model of Schaerer and de Koter (1997). 
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Figure 9: Calculated spectral energy distribution of a typical O star. The spectrum represents 
the result of the WM-basic model (Pauldrach etal. 2001), and the thick line represents the 
result of a corresponding COSTAR model (Schaerer and de Koter 1997). 

The spectral energy distributions have been calculated for two modes of cluster 
evolution, a continuous and an impulsive one. The different spectral evolutions of 
these modes are shown in Figure 10 with regard to the spectral range of 0.8-2.5 
Rydbergs. 

In the continuous mode, the star-formation rate is assumed to be constant with 

time - on time scales which are longer than the lifetimes of massive stars - and the 
cluster is assumed to form stars at a rate of 1 Mq per year. In the impulsive mode, 
all stars are formed "instantaneously", i. e., on time scales which are inuch shorter 
than the lifetimes of massive stars; the total mass of stars formed is lO'' M:v), and 
the cluster is assumed to evolve passively thereafter. Furthermore, the computations 
are based on a Salpeter initial mass function, where for each mode two cases are 
assumed, one with an upper mass limit of Mup = 120 Mq and the other one with an 
upper mass limit of A'/up = 30 Mq. 

Two striking effects are seen in Figure 10: 

1. In the continuous mode, the cluster ionizing spectrum becomes softer - i.e., 
steeper - between 1 and 10 Myr. This is due to the increase of the relative num- 
ber of late versus early-type OB stars during the evolution, since the less massive 
late-type stars have longer lifetimes. This effect is particularly noticeable for the 
cluster model with Mup = 120 Mq. 

2. In the impulsive mode, the magnitude of the Lyman break increases drastically 
because the massive stars disappear with time. 
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Figure 10: Spectral energy distributions of time-evolving stellar clusters. Two modes of 
cluster evolution are shown: continuous and impulsive. (Left): Evolution for continuous star 
formation (lasting 10^" yr) for two extreme stellar compositions. In the upper panel evolv- 
ing SEDs of a cluster with Mup = 120 Mq, and in the lower panel with Mup = 30 M© 
are shown. Note that at early phases the amount of both ionizing and non-ionizing photons 
increase due to the increasing number of new-formed stars. After a time period of ~ 5 Myr 
the number of ionizing photons remains roughly constant, because an equilibrium of stellar 
aging and stellar birth is achieved, whereas the number of non-ionizing photons remains pro- 
portional to age, because of accumulation effects of old stars which are the major contributors 
to this spectral range. (Right): Evolution of an impulsive star formation burst (lasting 10^ yr), 
for the same extreme stellar compositions as above. Note that the ionizing photons practically 
disappear at an age of ~ 10 Myr. Figure from Sternberg etal. 2002. 



The evolution of the photon emission rates for photons above the Lyman ((5h) 
and the He I (Que) ionization thresholds is shown in Figure 1 1 for both modes of the 
clusters. The key difference between the two modes leads to completely different 
shapes of the photon emission rates, which will be easy to distinguish in view of 
their photoionizing properties acting on their gaseous environments. It is also shown 
that increasing the upper IMF mass cut-off increases the photon emission rates, but 
does not change the shapes of the rates obtained during the evolution. The methods 
developed for determining extragalactic abundances and population histories from 
an analysis of H ii regions are thus very promising. 

From these investigations it is obvious that realistic spectral energy distributions 
of massive stars and stellar clusters are important for studies of their environments. 
The crucial question, however, whether the spectral energy distributions of massive 
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Figure 1 1 : Evolution of the Lyman and He I photon emission rates for a range of values of the 
upper IMF mass cut-off. The continuous star formation mode is considered on the left hand 
side, and the impulsive star formation burst on the right. Figure from Sternberg et al. 2002. 

stars are already realistic enough to be used for diagnostic issues of H ii regions has 
not been answered yet, since on the basis of the results obtained for the diagnostic 
diagram in Section 4. 1 it cannot be excluded that wrong fluxes could show the same 
improvements in the Ne++/Ne+ ratios just by chance. 

An answer to this question requires an ultimate test] 

This ultimate test is only provided by comparing the observed and synthetic UV 
spectra of the individual massive stars, based on the following reasons: 

1 . This test involves hundreds of spectral signatures of various ionization stages with 
different ionization thresholds which cover a large frequency range. 

2. Almost all of the ionization thresholds lie within the spectral range shortward of 
the Lyman ionization threshold (cf. Figure 12); thus, the ionization balances of all 
elements depend sensitively on the ionizing radiation throughout the entire wind. 

The ionization balance can be traced reliably through the strength and structure 
of the wind lines formed throughout the atmosphere. Hence, it is a natural and the 
only reliable step to test the quality of the ionizing fluxes by virtue of their direct 
product: the UV spectra ofO stars. 

But before we turn to this test we will continue our discussion with another as- 
trophysically important stellar object. 
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Figure 12: Spectral energy distribution shortward of the Lyman ionization threshold of an 
expanding NLTE model atmosphere for an O supergiant (Pauldrach etal. 2003). The small 
vertical bars indicate the ionization thresholds for all important ions; the ionization balance 
depends almost entirely on the ionizing flux, and this influence can be traced by the spectral 
lines in the observable part of the UV spectrum. 



5 Supemovae of Type la as Distance Indicators 

In order to complete the discussion of Hot Stars in the context of modern astronomy, 
we will concentrate now on the subject of Supernovae of Type la. With respect to 
diagnostic issues we will investigate the role of Supemovae of Type la as distance 
indicators. The context of this discussion regards, as a starting point, the current sur- 
prising result that distant SNe la at intermediate redshift appear fainter than standard 
candles in an empty Friedmann model of the universe. 

Type la supemovae, which are the result of the thermonuclear explosion of a 
compact low mass star, are currently the best known distance indicators. Due to their 
large luminosities (imax ~ 10^^ erg/s) they reach far beyond the local supercluster 
(cf. Saha etal. 2001 and references therein). It is thus not surprising that Type la 
supemovae have become the most important cosmological distance indicator over 
the last years, and this is not only due to their extreme brightness, but primarily due 
to their maximum luminosity which can be normalized by their light curve shape, 
so that these objects can be regarded as standard candles. This has been shown by 
observations of local SNe la which define the linear expansion of the local universe 
extremely well (Riess etal. 1999), which in tum is convincing proof of the accuracy 
of the measured distances. On the other hand, observations of distant SNe la (up to 
redshifts of about 1) have yielded strong evidence that the expansion rate has been 
accelerated 6 Gyr ago (cf. Riess etal. 1998 and Perlmutter etal. 1999). 
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Figure 13: Left: Hubble diagram of Type la Supemovae. In the upper panel the distance mod- 
ulus versus redshift is shown in the usual way, whereas the distance modulus has been normal- 
ized to an empty universe in the lower panels. Shown are the data of the High-z SN Search 
Team (filled squares; Riess etal. 1998) and the Supernova Cosmology Project (open squares; 
Perlmutter etal. 1999). In addition, lines of four cosmological models are drawn, where the 
full line corresponds to a cosmological model for an empty universe and the dotted line to a 
flat universe. Figure from Leibundgut 2001. 

Right: Cosmological diagram showing the likelihood region as defined by SNe la in the 
versus S^Mass plane. The results of 79 SNe la have been included - 27 local and 52 distant 
ones ~ from the sources above. Contours indicate 68 %, 90 %, 95 %, and 99 % probability, 
and the line for a flat universe is indicated. Note that the SN la luminosity distances clearly 
indicate an accelerated expansion of the universe. Figure from Perlmutter etal. 1999. 



This surprising result becomes evident if we look at the Hubble diagram in the 
form of distance modulus versus redshift (cf. Figure 13). As is shown in the diagrams 
of the lower panels, which are normalized to a cosmological model for an empty 
universe, most SNe la at intermediate redshift are positioned at positive values of 
the normalized distance modulus; thus, the distant supernovae appear fainter than 
what would be expected in a empty universe. This means that deceleration from 
gravitational action of the matter content does not take place. Moreover, the SNe la 
appear even more distant indicating an accelerated expansion over the last 6 Gyr. 

The implications of this are immediately recognized in the cosmological diagram 
on the right-hand side of Figure 13, where n\ versus f^Mass is plotted. The SN la 
results obviously exclude any world model with matter but without a a cosmological 
constant. Hence, the data show the need for an energy contribution by the vacuum. 
Thus, the SN luminosity distances indicate accelerated expansion of the universe! 

However, this is not the only interpretation of the result obtained. There are other 
astrophysical explanations, such as obscuration by intergalactic dust or evolution of 
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Figure 14: Observed HST spectrum of a standard Supernova of Type la at early phase - 
SN 1992A from 24 Jan. 1992. Figure from Kirshner et al. 1993. 



supernovae of Type la (cf. Leibundgut 2001). The latter point refers primarily to 
the absolute peak luminosity of SNe la which might have changed between the local 
epoch and a redshift of 1.0, especially in the case of a decrease: then the fainter 
supernovae would not be a signature of larger distances, but rather of an evolution 
of these objects. Although present SN la models do not favor such an evolution, it 
cannot be excluded, because neither the explosions nor the radiation emerging from 
the atmosphere are understood in detail yet. Consequently, serious caveats for the 
cosmological interpretation of distant supernovae exist. 

Thus, we are faced with the question: Are SNe la standard candles independent 
of age, or is there some evolution of the SN luminosity with age? Among other 
things spectroscopy is certainly a powerful tool to obtain an answer to this question 
by searching for spectral differences between local and distant SNe la. 

Unfortunately, at present, we do not have a clear picture of the exact physical 
processes which take place in a SN la explosion and how the radiation released from 
it should be treated (for a recent review see Hillebrandt and Niemeyer 2000). In 
order to get an impression of the basic physics involved that affect the atmospheric 
structure of a SN la, a quick glance at an observed spectrum will help to point out 
the necessary steps for their further analysis. 

Figure 14 shows a typical SN la spectrum at early epochs (< 2 weeks after 
maximum). The most striking features of this spectrum are the characteristic P Cygni 
line profiles which are quite similar to the signatures of O stars. But compared to the 
latter objects there are also important differences: the broad lines indicate that the 
velocities of the SN la ejecta are almost a factor of 10 larger (up to 30000 km/s) 
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and, as a second important point, SN la spectra contain no H and He lines. Instead, 
prominent absorption features of mainly intermediate-mass elements (Si ii, O i, S ii, 
Ca II, Mg II, . . . ) embedded in a non-thermal pseudo-continuum iire observed. 

To answer the question of whether SN la are standard candles in a cosmological 
sense, realistic models and synthetic spectra of Type la Supemovae are obviously 
required for the analysis of the observed spectra. As just shown, however, these 
models will necessarily need to be based on a similar serious approach for expanding 
atmospheres, characterized by non-equilibrium thermodynamics, as is the case for 
Hot Stars in general. 

It will be shown in Section 7.1 that such spectra are in principle already available. 

6 Concept for Consistent Models 
of Hot Star Atmospheres 

In order to determine stellar abundances, parameters, and physical properties (and 
from these, obtain realistic spectral energy distributions) of Hot Stars via quantita- 
tive UV spectroscopy, a principal difficulty needs to be overcome: the diagnostic 
tools and techniques must be provided. This requires the construction of detailed 
atmospheric models. 

As has been demonstrated in the previous sections, such a tool has not been 
made obsolete by general development in astrophysics. On the contrary, it is becom- 
ing more and more relevant to current astronomical research as spectral analysis of 
hot luminous stars is of growing astrophysical interest. Thus, continuing effort is ex- 
pended to develop a standard code in order to provide the necessary diagnostic tool. 
In the following part of this section we describe the status of our continuing work to 
construct realistic models for expanding atmospheres. 

Before we focus on the theory in its present stage, we should mention previous 
fundamental work which turned out to be essential in elaborating the theory. In this 
context I want to emphasize publications which refer to key aspects of theoretical 
activity. The starting point of the development of the radiation-driven wind theory 
is rooted in a paper by Milne (1926) more than 70 years ago. Milne was the first 
to realize that radiation could be coupled to ions and that this process subsequently 
may eject the ions from the stellar surface. The next fundamental step goes back 
to Sobolev (1957), who developed the basic ideas of radiative transfer in expanding 
atmospheres. 

Radiation pressure as a driving mechanism for stellar outflow was rediscovered 
by Lucy and Solomon (1970) who developed the basis of the theory and the first 
attempt to its solution. The pioneering step in the formulation of the theory in a 
quasi self-consistent manner was performed by Castor, Abbott, and Klein (1975). 
Although these approaches were only qualitative (due to many simplifications), the 
theory was developed further owing to the promising results obtained by these au- 
thors. Regarding key aspects of the solution of the radiative transfer the work of Ry- 
bicki (1971) and Hummer and Rybicki (1985) has to be emphasized. With respect to 
the soft X-ray emission of O stars (detected by Seward et al. 1979 and Harnden etal. 
1979), Cassinelli and Olson (1979) have investigated the possible influence of X-rays 
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Figure 15: Schematic sketch of the basic equations and the non-linear system of integro- 
diffcrcntial equations that form the basis of stationary models of hot star atmospheres - for 
details see text. Figure from Pauldrach etal. 1994b). 
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on the ionization structure. Finally, Lucy and White (1980) and Owocki, Castor, and 
Rybicki (1988) have distinguished themselves with basic theoretical investigations 
of time-dependent radiation hydrodynamics which describe the creation and devel- 
opment of shocks (for more details about the role of X-rays in the atmospheres of 
Hot Stars see the reviews of Pauldrach etal. 1994b and Kudritzki and Puis 2000). 

6.1 The general method 

The basis of our approach in constructing detailed atmospheric models for Hot Stars 
is the concept of homogeneous, stationary, and spherically symmetric radiation 
driven winds, where the expansion of the atmosphere is due to scattering and ab- 
sorption of Doppler-shifted metal lines. Although these approximations seem to be 
quite restrictive, it has already been shown that the time-averaged mean of the ob- 
served UV spectral features can be described correctly by such a method (Pauldrach 
etal. 1994, 1994b). 

The primary output of this kind of model calculation are the spectral energy 
distributions and synthetic spectra emitted by the atmospheres of hot stars. As these 
spectra consist of hundreds of strong and also weak wind-contaminated spectral lines 
which form the basis of a quantitative analysis, and as the spectral energy distribution 
from hot stars is also used as input for the analysis of emission line spectra which 
depend sensitively on the structure of the emergent stellar flux (cf. Section 4), a so- 
phisticated and well-tested method is required to produce these data sets accurately. 

However, developing such a method is not straightforward, since modelling the 
atmospheres of hot star involves the replication of a tightly interwoven mesh of phys- 
ical processes: the equations of radiation hydrodynamics including the energy equa- 
tion, the statistical equilibrium for all important ions with detailed atomic physics, 
and the radiative transfer equation at all transition frequencies have to be solved si- 
multaneously. Figure 15 gives an overview of the physics to be treated. 

The principal features are: 

(1) The stellar parameters Teff (effective temperature), log 5 (logarithm of photo- 
spheric gravitational acceleration), ii* (photospheric radius defined at a Rosse- 
land optical depth of 2/3) and Z (abundances) have to be pre-specified. 

(2) The hydrodynamic equations are solved (r is the radial coordinate, g the mass 
density, v the velocity, p the gas pressure and M the mass loss rate). The 
crucial term is the radiative acceleration ^rad with minor contributions from 
continuous absorption and major contributions from scattering and line absorp- 
tion. For each line the oscillator strengths /;„, the statistical weights gi, Qu, and 
the occupation numbers nu n„ of the lower and upper level enter the equations, 
together with the frequency and angle integral over the specific intensity 7^ and 
the line broadening function ipi, accounting for the Doppler effect. 

(3) The occupation numbers are determined by the rate equations containing colU- 
sional (Cy ) and radiative {Rij) transition rates. Low-temperature dielectronic 
recombination is included and Auger ionization due to K-shell absorption (es- 
sential for C, N, O, Ne, Mg, Si, and S) of soft X-ray radiation arising from 
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shock-heated matter is taken into account using detailed atomic models for all 
important ions. Note that the hydrodynamical equations are coupled directly 
with the rate equations. The velocity field enters into the radiative rates while 
the density is important for the colhsional rates and the equation of particle 
conservation. On the other hand, the occupation numbers are crucial for the 
hydrodynamics since the radiative line acceleration dominates the equation of 
motion. 

(4) The spherical transfer equation which yields the radiation field at every depth 
point, including the thermalized layers where the diffusion approximation is 
applied, is solved for the total opacities {Kv) and source functions {S^,) of all 
important ions. Hence, the influence of the spectral lines - the strong EUV 
line blocking - which affects the ionizing flux that determines the ionization 
and excitation of levels, is taken into account. Note that the radiation field is 
coupled with the hydrodynamics (grad) and the rate equations (Rij). 

Moreover, the shock source functions (S^) produced by radiative coohng zones 
which originate from a simulation of shock heated matter, together with K- 
sheU absorption (k^), are also included in the radiative transfer. The shock 
source function is incorporated on the basis of an approximate calculation of 
the volume emission coefficient (A^) of the X-ray plasma in dependence of the 
velocity-dependent post-shock temperatures and the filling factor (/). 

(5) The temperature structure is determined by the microscopic energy equation 
which, in principle, states that the luminosity must be conserved. Line blan- 
keting effects which reflect the influence of line blocking on the temperature 
structure are taken into account. 

(6) The iterative solution of the total system of equations yields the hydrodynamic 

structure of the wind (i. e., the mass-loss rate M and the velocity structure v{r)) 
together with synthetic spectra and ionizing fluxes. 

Essential steps which form the basis of the theoretical framework developed are de- 
scribed in Pauldrach etal. 1986, Pauldrach 1987, Pauldrach and Herrero 1988, Puis 
and Pauldrach 1990, Pauldrach et al. 1994, Feldmeier et al. 1997, and Pauldrach et al. 
1998. 

The effect complicating the system the most is the overlap of the spectral lines. 
This effect is induced by the velocity field of the expanding atmosphere which shifts 
at different radii up to 1000 spectral fines of different ions into the line of sight 
at each observer's frequency. Since the behavior of most of the UV spectral Unes 
additionally depends critically on a detailed and consistent description of the corre- 
sponding effects of line blocking and line blanketing (cf. Pauldrach 1987), special 
attention has to be given to the correct treatment of the Doppler- shifted line radiation 
transport of all metal Unes in the entire sub- and supersonically expanding atmo- 
sphere, the corresponding coupling with the radiative rates in the rate equations, and 
the energy conservation. Another important point to emphasize concerns the atomic 
data, since it is quite obvious that the quality of the calculated spectra, the multilevel 
NLTE treatment of the metal ions (from C to Zn), and the adequate representation 
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of line blocking and the radiative line acceleration depends crucially on the quaUty 
of the atomic models. Thus, the data have to be continuously improved whenever 
significant progress of atomic data modelling is made. Together with a revised inclu- 
sion of EUV and X-ray radiation produced by cooling zones which originate from 
the simulation of shock heated matter, these improvements have been implemented 
and described in a recent paper (Pauldrach et al. 2001). Very recently we have addi- 
tionally focused on a remaining aspect regarding hydrodynamical calculations which 
are consistent with the radiation field obtained from the line-overlap computations, 
and we have incorporated this improvement into our procedure (see Pauldrach and 
Hoffmann 2003). 

This solution method is in its present stage already regarded as a standard proce- 
dure towards a realistic description of stationary wind models. Thus, together with 
an easy-to-use interface and an installation wizard, the program package WM-basic 
has been made available to the community. The package can be downloaded from 
the URL given on the first page. 

In the next section we will investigate thoroughly whether these kind of models 
already lead to reaUstic results. 

7 Synthetic Spectra and Models 
of Hot Star Atmospheres 

The first five sections of this review inevitably lead to the realization that line fea- 
tures of expanding atmospheres are one of the most important astronomical manifes- 
tations, since they can easily be identified even in spectra of medium resolution in 
individual objects like supernovae or in integrated spectra of starburst regions even at 
significant redshift. As a consequence, they can be used to provide important infor- 
mation about the chemical composition, the extragalactic distance scale, and several 
other properties of miscellaneous stellar populations. All that is required are hydro- 
dynamic NLTE model atmospheres that incorporate the effects of spherical extension 
and expansion as realistic and consistent as possible. After almost two decades of 
work these model atmospheres now seem to be available, and in the following we 
will investigate thoroughly whether basic steps towards realistic synthetic spectra 
for Supernovae of Type la can already be presented, and whether UV diagnostic 
techniques applied to two of the most luminous O supergiants in the Galaxy lead to 
synthetic spectra which can be regarded as a measure of excellent quality. 

7.1 Synthetic and Observed Spectra 
of Supernovae of Type la 

Prominent features in the spectra of Type la supernovae in early phases are the char- 
acteristic absorption lines resulting from low ionized intermediate-mass elements, 
such as Sin, Oi, Can, Mgii, cf. Figure 17. These absorption lines show char- 
acteristic line shapes due to Doppler-broadening resulting from the large velocity 
gradients. The formation of these lines results in a pseudo-continuum that is set up 
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Figure 16: Synthetic NLTE spectra of a Super- 
nova of Type la. Top: The dotted hnes show the 
spectrum from a model without line blocking 
(left panel) and from one where line blocking 
is treated consistently (right panel), both com- 
pared to the observed spectrum of SN 1992 A 
(solid line). (Figure from Pauldrach et al. 1996.) 
Left: The two curves in the lower part show 
a synthetic spectrum for a W7 model at 23 
days compared to the observed spectrum of 
SN 1992A. (Figure from Nugent etal. 1997.) 
The observed spectrum is from Kirshner etal. 
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by the enormous number of these lines. Supernovae appear in this 'photospheric' 
epoch for about one month after the explosion and the spectra at this epoch con- 
tain useful information on the energetics of the explosion (luminosity, velocity- and 
density-structure) and on the nucleosynthesis in the intermediate and outer part of the 
progenitor star. Therefore, the primary objective in order to analyze the spectra is to 
construct consistent models which link the results of the nucleosynthesis and hydro- 
dynamics obtained from state-of-the-art explosion models (cf. Reinecke etal. 1999) 
with the calculations of light curves and synthetic spectra of SN la. In a first step, 
however, this is done by using the standard hydrodynamical model W7 by Nomoto 
etal. 1984. Figure 16 exemplifies the status quo of the synthetic Spectra of Type la 
Supernovae at early phases obtained from this first step, compared to observations. 

Obviously, the first comparison shown (top left panel of Figure 16) does not at 
all represent a reaUstic atmospheric model. The synthetic spectrum shown has been 
calculated for a model without line blocking, i. e., only continuum opacity has been 
considered in the NLTE calculation. Thus, the drastic effects of line blocking on 
the ionization and excitation and the emergent flux can be verified by comparing the 
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Figure 17: Comparison of the most recent synthetic NLTE spectrum at early phases and the 
observed spectrum of SN 1992A, illustrating that the method used is already on a quantitative 
level. The rest wavelengths of various metal lines are indicated by vertical bars. Figure from 
Sauer and Pauldrach 2002. 

calculated flux with the observed flux, where the huge difference in the UV is due 
to the missing contribution of the lines to the opacity. It is an essential by-product 
of this result that the observed line features in the optical are also not reproduced: 
The two strongest features in the synthetic spectrum are due to Fe iii lines, indicating 
that the ionization is too high because of the excessive UV flux. This behavior nicely 
illustrates that the ionization balance depends almost entirely on the ionizing UV 
flux and that this influence can be traced by the spectral lines in the observable part 
of the spectrum, as has been stated in Section 4.2. 

The fact that the difference in the UV observed for this model is indeed due to 
line blocking becomes obvious if this important effect is properly taken into account, 
as in the case of the model shown in the upper right panel of Figure 16. The agree- 
ment of the resulting spectrum of this model with the observations is quite good, 
demonstrating that line blocking alone can diminish the discrepancy between the 
synthetic and the observed spectra. Both the flux level and most line features are 
now reasonably well reproduced throughout the spectrum (cf. Pauldrach et al. 1996). 
Additionally, in the lower left panel of Figure 16 a spectral fit to the same super- 
nova observation by Nugent et al. (1997) is shown, which is based on a completely 
independent approach to the theory. The resulting spectrum is also quite reasonable. 

Finally, Figure 17 shows a recently calculated synthetic spectrum that is based 
on the important improvements of the theory described in the previous section, com- 
pared to the observed "standard" supernova spectrum (SN 1992A). The synthetic 
spectrum reproduces the observed spectrum quite well, and the overall impression of 
this comparison is that the method used is already on a quantitative level, indicating 
that the basic physics is treated properly (cf. Sauer and Pauldrach 2002). 
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Although there is still need for improvements regarding, for instance, a proper 
treatment of the 7 energy deposition in the outermost layers and the use of current 
hydrodynamic explosion models, these models form the basis for the primary objec- 
tive of the subject, i. e., to search for spectral differences between local and distant 
SNe la. Concerning the faciUties of spectral analysis we are thus close to the point 
where we can tackle the question of whether SNe la are standard candles in a cos- 
mological sense. 

7.2 Detailed Analyses of Massive O Stars 

The objectives of a detailed comparison of synthetic and observed UV spectra of 
massive O stars are manifold. Primarily it has to be verified that the higher level 
of consistent description of the theoretical concept of line blocking and blanketing 
effects and the involved modifications to the models leads to changes in the Une 
spectra with much better agreement to the observed spectra than the previous, less 
elaborated and less consistent models. Secondly it has to be shown that the stellar 
parameters, the wind parameters, and the abundances can be determined diagnosti- 
cally via a comparison of observed and calculated high-resolution spectra covering 
the observable UV region. And finally, the quality of the spectral energy distributions 
has to be verified. The latter point, however, is a direct by-product of the other objec- 
tives if the complete observed high-resolution UV spectra are accurately reproduced 
by the synthetic ones, as illustrated in Section 4.2 

In the following, the potential of the improved method will be demonstrated by 
an application to the 04 1(f) star C, Puppis and the 09.5 la star a Cam. 

7.2.1 UV Analysis of the Hot O Supergiant C Puppis 

Referring to the beginning of the review, we will now start to examine carefully the 
old-fashioned UV spectrum of C, Puppis. To put this into perspective, we will also 
briefly review the improvements of UV fine fits foUowing from gradual improve- 
ments of the methods used. 

The first serious attempt to analyze the UV spectrum of this standard object goes 
back to Lamers and Morton in 1976. From a present point of view their model can 
not be regarded as a sophisticated one, because they assumed the complete model 
structure - the dynamical structure, the occupation numbers, and the temperature 
structure - and the fine radiative transfer was treated in an approximative way. Nev- 
ertheless, we learned from this kind of spectrum synthesis that a solution to the 
problem is feasible, and that at least the dynamical parameters - the mass-loss rate 
(M) and the terminal velocity {Voo) - can, in principle, be determined from UV lines 
quite accurately. The UV line fits of Hamann in 1980 have been based on the same 
assumptions apart from a very detailed treatment of the radiative line transfer. As 
a result of this the calculated resonance lines of some ions were already quite well 
in agreement with the observed ones. But the real conclusion to be drawn from 
these comparisons is that a detailed analysis requires a more consistent treatment of 
expanding atmospheres. 



Hot Stars: Old-Fashioned or Trendy? 



11 




Figure 18: Calculated and observed UV spectrum of <^ Puppis. The observed spectrum shows 
the Copernicus and lUE high-resolution observations, and the calculated spectrum belongs to 
the final model of Pauldrach etal. 1994. 



The situation improved considerably with the first attempts to find a consistent 
solution, meaning that the equations of hydrodynamics, non-equilibrium thermody- 
namics, and radiative transfer have been solved in a consistent way. Pauldrach (1987) 
introduced for the first time a full NLTE treatment of the metal lines driving the 
wind, and Puis (1987) investigated the important effect of multiple photon momen- 
tum transfer through line overlaps caused by the velocity-induced Doppler- shifts for 
applications in stellar wind dynamics. From this procedure, dynamical parameters, 
constraints on the stellar parameters, and, as has been verified by a comparison of a 
sample of calculated spectral lines with the observed spectrum, an ionization equilib- 
rium and occupation numbers which were close to a correct description have already 
been obtained. But, among other approximations, the models stiU suffered from the 
neglect of radiation emitted from shock cooling zones and from a very approximate 
treatment of hne blocking. 

The further steps to reproduce most of the important observed individual line 
features in the U V required a lot of effort in atomic physics, in improved NLTE mul- 
tilevel radiative transfer, and in spectrum synthesis techniques. An important step 
towards this objective was the paper by Pauldrach etal. (1994), where for the first 
time stellar parameters and abundances have been determined from individual UV 
hne features. Figure 18 shows the corresponding synthetic spectrum that was al- 
ready in overall agreement with the observations. Nevertheless, the treatment was 
still affected by a number of severe approximations regarding especially the neglect 
of the hne blanketing effect, an approximate treatment of the line radiative transfer. 
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Figure 19: Calculated and observed UV spectrum for Puppis. The observed spectrum 
shows the Copernicus and IUE high-resolution observations, and the calculated state-of-the- 
art spectra represent the final models of Pauldrach et al. 2003. For a discussion see text. 
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Figure 20: Calculated state-of-the-art spectral energy distribution of the final model of Pup- 
pis of Pauldrach etal. 2003. Left: Linear scale for the flux. Right: Logarithmic scale for the 
flux. 

and a simple treatment of emitted radiation from shock cooling zones. Thus, we 
conclude that a more consistent treatment of expanding atmospheres was still needed 
for a detailed analysis. 

This more consistent treatment of expanding atmospheres that has now been ap- 
plied and fixes the status quo is based on the improvements discussed in Section 6.1 
(for a more comprehensive discussion see Pauldrach etal. 2001 and Pauldrach and 
Hoffmann 2003). As is shown in Figure 19, the calculated synthetic spectrum is quite 
well in agreement with the spectra observed by lUE and Copernicus. The small dif- 
ferences observed in the Si iv and the N iv lines just reflect a sensitive dependence 
on the parameters used to describe the shock distribution. This is verified by a com- 
parison of the two panels of Figure 19, where just the shock-distribution has slightly 
been changed within the range of uncertainty of the corresponding parameters. Not 
only have the stellar and wind parameters of this object been confirmed by the model 
on which the synthetic spectrum is based, but also the abundances of C, N, O, P, Si, 
S, Fe, and Ni have been determined. We thus conclude that the present method of 
quantitative spectral UV analysis of Hot Stars leads to models which can be regarded 
as being realistic. 

Consequently, we consider these kind of quantitative spectral UV analyses as the 
ultimate test for the accuracy and the quality of theoretical ionizing fluxes (cf . Fig- 
ure 20), which can thus be used as spectral energy distributions for the analysis of 
H II regions. 

7.2.2 UV Analysis of the Cool O Supergiant a Cam 

In Figure 21 it can be seen that the method also works for cool O supergiants Uke 
a Cam. We recognize again that the synthetic spectrum is quite well in agreement 
with the rUE and Copernicus spectra (for a discussion see Pauldrach etal. 2001). 
In the uppermost panel of the spectrum, simulations of interstellar lines have addi- 
tionaUy been merged with the synthetic spectrum in order to disentangle wind lines 
from interstellar lines (for a discussion of this point see Hoffmann 2002). In this case 
also, our new realistic models allowed to determine stellar parameters, wind param- 
eters, and abundances from the UV spectra alone. (The operational procedure of the 



30 



A. W.A. Pauldrach 




1500 1550 1600 1650 1700 1750 1800 1850 1900 

Wavelength (A) 

Figure 21: Comparison of the best model of Pauldrach etal. 2001 with spectra of a Cam 
observed by lUE and Copernicus, demonstrating the quality that can be achieved with the new 

model generation. 

method begins with realistic estimates of R, Toft, M, and a set of abundances. With 
these, the model atmosphere is solved and the velocity field, the mass loss rate M, 
and the synthetic spectrum is calculated. The parameters are adjusted and the process 
is repeated until a good fit to all features in the observed UV spectrum is obtained. 
For a compact description of this procedure see Pauldrach etal. 2002). It turned out 
that the effective temperature can be determined to within a range of ±1000 K and 
the abundances to at least within a factor of 2. 

7.3 Wind properties of massive O stars 

As a last point we want to illustrate the significance of the dynamical parameters of 
radiation-driven winds. The intrinsic significance is quite obvious: it is the consistent 
hydrodynamics which provides the link between the stellar parameters (Teff , M, R) 
and the wind parameters (voo, M). Thus, the appearance of the UV spectrum is 
determined by the interplay of the NLTE model and the hydrodynamics. (As was 
discussed in Section 6.1, the hydrodynamics is controlled by the line force, which 
is primarily determined by the occupation numbers, and the radiative transfer of the 
NLTE model, but the hydrodynamics in turn affects the NLTE model via the density 
and velocity structure.) 

A tool for illustrating the significance of the dynamical parameters is offered 
by the so called wind-momentum-luminosity relation. This relation is based on two 
important facts: The first one is that, due to the driving mechanism of hot stars. 
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Figure 22; The wind-momentum-luminosity relation for massive O stars and Central Stars 
of Planetary Nebulae (CSPNs). Circles designate the O star analysis based on Ha profiles by 
Puis etal. 1996, and squares that of CSPNs by Kudritzki etal. 1997. Figure from Kudritzki 
and Puis 2000. 



the mechanical momentum of the wind flow (voqM) is mostly a function of photon 
momentum (L / c) and is therefore related to the luminosity. The second one is that 
the expression VaoMR^f^ is an almost directly observable quantity. As the shapes 
of the spectral lines are characterized by the strength and the velocity of the outflow, 
the first term of this expression, the observed terminal velocity, can be measured 
directly from the width of the absorption part of the saturated UV resonance lines. 
Deducing the mass-loss rate from the Une profiles is, however, more complex, since 
this requires the calculation of the ionization balance in advance. On the other hand, 
the advantage is that optical lines like Ha can be used for this purpose. Thus, the 
product of the last two terms of the expression VaoMR^^"^ directly follows from a 
line fit of Ha (cf. Puis etal. 1996). The dynamical parameters obtained in this way 
are usually designated as observed wind parameters. Figure 22 shows that the wind- 
momentum-luminosity relation indeed exists for massive O stars, as they follow 
the linear relation predicted by the theory, in a first approximation and for fixed 
metallicities. With regard to the spread in wind momenta found by Kudritzki et al. 
(1997) for the Central Stars of Planetary Nebulae (lower part of Figure 22), Pauldrach 
etal. (2002, 2003) have given a solution which, however, is unlikely to be believed 
yet by the part of the community working on stellar evolution. 

Moreover, it is important to note that as a by-product of this relation it can be 
used as an independent tool for measuring extragalactic distances up to Virgo and 
Fornax, since it is independent of the observationally unknown masses (for a recent 
and comprehensive review on this subject see Kudritzki and Puis 2000). 

Figure 23 shows that the observed behavior of the wind momentum luminosity 
relation is represented quite well by our improved realistic models, particularly in 
the case of the two supergiants (C Puppis and a Cam) analyzed in Section 7.2. But, 
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Figure 23: The wind-momentum-luminosity relation for massive O stars and CSPNs as in 
Figure 22. Also plotted as filled squares are the calculated wind momenta for five massive 
O stars; filled circles designate the corresponding observed values (cf. Pauldrach etal. 2002, 
2003). 



as already pointed out, this relation is independent of the stellar mass. Thus, in 
order to verify the statement of the previous section that the observed stellar and 
wind parameters are confirmed by the present models, we also have to investigate 
the relations of the individual dynamical parameters. 

For this investigation we need to use as input for our models the same stellar 
parameters as have been used to obtain the observed wind parameters. On the basis 
of this requirement. Figure 24 (upper panel) shows that the observed and predicted 
values of the terminal velocities (Voo) are in agreement to within 10 %. Since Voo is 
proportional to the escape velocity (Uesc) 

/2GM y/' 

which strongly depends on the mass of the objects, the mass is determined very 
accurately by the predicted values (G is the gravitational constant and F is the ratio 
of radiative Thomson acceleration to gravitational acceleration). For the mass-loss 
rates we found agreement to within a factor of two, as shown in Figure 24 (lower 
panel). Due to the strong relation between the mass-loss rates and the luminosities 
(M OC L), the luminosities can thus be precisely determined. 

Therefore, computing the wind dynamics consistently with the NLTE model per- 
mits not only the determination of the wind parameters from given stellar parameters, 
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Figure 24: Top: Terminal velocities as a function of effective temperature for massive 0-stars 
of the same sample as in Figure 23. Bottom: Mass-loss rates as a function of effective temper- 
ature for the same sample. The calculated values for the five massive O stars are designated 
by filled squares; filled circles designate the corresponding observed values (cf. Pauldrach and 
Hoffmann 2003 concerning C, Puppis; the other objects have been treated by Hoffmann and 
Pauldrach 2002). 
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but conversely makes it possible to obtain the stellar parameters from the observed 
UV spectrum alone. As stated above, this means that, in principle, the stellar param- 
eters can immediately he read off by simply comparing an observed UV spectrum to 
a proper synthetic spectrum. However, the whole procedure is not always an easy 
task, since in most cases a comprehensive grid of models and some experience will 
be required to do so. Although the idea itself is not new (cf. Pauldrach et al. 1988 and 
Kudritzki etal. 1992), only the new generation of models has reached the degree of 
sophistication that makes such a procedure practicable instead of purely academic. 
The corresponding conclusion is that realistic models are characterized by at least a 
quantitative spectral UV analysis calculated together with a consistent dynamics. 



8 Conclusions and Outlook 

The need for detailed atmospheric models of Hot Stars has been motivated in 
depth and a diagnostic tool with great astrophysical potential has been presented. 

It has been shown that the models of the new generation are realistic. They are 
realistic with regard to a quantitative spectral UV analysis calculated, for the case 
of O stars, along with consistent dynamics, which, in principle, allows to read off 
the stellar parameters by comparing an observed UV spectrum to a suitable synthetic 
spectrum. The new generation of models has reached a degree of sophistication that 
makes such a procedure practicable. 

The astronomical perspectives are enormous. Of course, there is still a large 
amount of hard and careful astronomical work that needs to be done. The diagnostic 
tool for expanding model atmospheres is in our hands, but this is (as usual) just the 
starting point; further elaboration, refinements, and modifications are required before 
the results are quantitatively completely reliable. 

But this is not the handicap for the future. The present handicap are future ob- 
servations. What is most urgently needed is the Next Generation Space Telescope. 

With respect to the topics discussed in this review this telescope is needed for a 
variety of reasons: 

To tackle the question whether Supernovae la are standard candles in a cosmo- 
logical sense realistic models and synthetic spectra of Type la Supernovae require 
spectral observations of objects at mid-redshift and sufficient resolution. The con- 
text of this question is the current surprising result that distant SNe la at intermediate 
redshift appear fainter than standard candles in an empty Friedmann model. Con- 
sequently, the current SN-luminosity distances indicate an accelerated expansion of 
the universe. 

In order to be able to make quantitative predictions about the influence of very 
massive, extremely metal-poor Population HI stars on their galactic and intergalac- 
tic environment one primarily needs observations which can be compared to the pre- 
dicted flux spectra that are already available for zero metalhcity and can, in principle, 
be produced for metallicities different from zero. Observations with the Next Gen- 
eration Space Telescope (NGST) of distant stellar populations at high redshifts will 
give us the opportunity to deduce the primordial IMF, thus allowing a quantitative 
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investigation of the ionization efficiency of a Top-heavy IMF via realistic spectral 
energy distributions of these very massive stars. 

NGST observations are also important for determining extragalactic abundances 
and population histories of starburst galaxies from an analysis of H ii regions. This 
task also requires energy distributions of time-evolving stellar clusters, which are cal- 
culated on the basis of a grid of spectral energy distributions of O and early B-type 
stars. A crucial point with respect to this is whether the spectral energy distribu- 
tions of massive stars are already reaUstic enough to be used for diagnostic issues of 
H II regions. It has been shown in this review that the ultimate test is provided by a 
comparison of observed and synthetic UV spectra of individual massive stars, since 
the ionization balance can be traced reliably through the strength and structure of the 
wind lines formed throughout the atmosphere. 

Furthermore, NGST is important for directly exploiting the diagnostic perspec- 
tives of galaxies with pronounced current star formation. It has been demonstrated 
that massive stars dominate the UV wavelength range in star-forming galaxies, and 
that therefore the UV-spectral features of massive O stars can be used as tracers of 
age and chemical composition of starburst galaxies even at high redshift. This is in 
particular the case when the flux from these galaxies is amplified by gravitational 
lensing through foreground galaxy clusters; corresponding observations render the 
possibility of determining metaUicities of starbursting galaxies in the early universe 
via realistic UV spectra (in the rest frame) of massive O stars. 

As a final remark it is noted that the solution method of stationary models for ex- 
panding atmospheres is in its present stage already regarded as a standard procedure 
towards a realistic description. Thus, together with an easy-to-use interface and an 
installation wizard, the program package WM-basic has been made available to the 
community and can be downloaded from the author's home page. 
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